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1—5 ; 1 Introduction 

7—1 . En réponse à la demande croissante en une densité d'informations stockées et traitées de plus en plus élevée, 

r^i, un effort important est consacré à la conception de nouveaux matériaux pour le stockage et le traitement 
D ! de l'information à l'échelle moléculaire. Les systèmes de métaux de transition occupent dans ce domaine 

'5 | une place essentielle en raison de leurs propriétés photophysiques et photo chimiques très variées, et de la 
^ | possibilité de coupler celles-ci à leurs propriétés magnétiques. C'est précisément le cas des complexes à fran- 
co ' sition de spin et leurs composés qu'il est possible de faire évoluer réversiblement entre leur état électronique 
S l fondamental bas-spin et leur état électronique excité métastable haut-spin, de propriétés magnétiques et 

'"O ' optiques distincts. Le changement d'états de spin peut être induit par variation de la température, par 
q ■ application d'une pression externe ou d'un champ magnétique, et par excitation lumineuse. Du fait de la 

i vaste étendue de fonctionnalités auxquelles ces systèmes photomagnétiques bistables permettraient d'accé- 
. der, ceux-ci sont considérés comme les archétypes des commutateurs optomagnétiques. Ils sont ainsi l'objet 

r ^ , de nombreuses études multidisciplinaires [T]. 

Nos travaux sur les complexes à transition de spin ou "apparentés" s'inscrivent dans le cadre de l'étude 
| de la commutation photoinduite. La compréhension approfondie de ce phénomène requiert de pouvoir 

CN | répondre aux questions clés suivantes : 

7— I ■ 

! (Qi) Quelles sont les caractéristiques du système considéré dans les états entre lesquels on souhaite le faire 
\ évoluer? (structures, énergies, durées de vie, etc. . .) 

• i-H ', (Q2) Comment le système interagit-il avec la lumière? (énergies d'excitation, et moments de transition 
^ \ dipolaires et/ou magnétiques) 

S; 

(Q3) Comment évolue- t-il après la photo-excitation? 

Toutefois, apporter des réponses précises à ces trois questions reste un défi majeur que nous essayons de 
contribuer à relever par nos travaux basés sur la DFT (pour l'anglais density-functional theory) et la TDDFT 
(pour l'anglais time- dépendent DFT). Dans ce chapitre, nous donnons un aperçu des problématiques liées 
à l'application de la DFT et de la TDDFT à l'étude des complexes de métaux de transition, avec des 
exemples principalement tirés de nos travaux, et, pour plus de détails, renvoyons le lecteur aux références 
données. L'organisation du chapitre est la suivante. Dans la prochaine section, nous faisons une brève revue 
méthodologique, tout en mettant l'emphase sur la DFT et la TDDFT. Dans la Sec. [3l nous faisons le point 
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sur la performance des différentes fonctionnelles de la DFT et montrons que les fonctionnelles modernes 
rivalisent en qualité avec les méthodes ab initio mais avec un coût computationnel moindre. Dans la Sec. [4j 
nous considérons l'apport de la TDDFT à l'étude des états excités des complexes. Dans la Sec. [5j nous 
abordons le sujet important des effets d'environnement sur le comportement des complexes à transition de 
spin. Finalement, dans la Sec. |6j nous concluons et essayons d'anticiper l'évolution future du domaine. 



2 Méthodologie : Trois points de Pauling 

" Essentially, ail models are wrong, but some are useful." [2] Cette phrase du statisticien américain George 
E. P. Box, prise hors contexte, nous rappelle que nos modèles restent toujours des modèles et par conséquent, 
qu'en dépit de nos meilleurs efforts, ils ont tous des défauts allant de pair avec leurs avantages. Bien souvent, 
lors du développement d'un modèle, la qualité de l'accord théorie-expérience n'augmente pas de manière 
monotone avec le degré d'élaboration du modèle mais a plûtot tendance à osciller avec, à certains niveaux 
de sophistication, des minima sur la courbe « erreur théorie-expérience » / « sophistication du modèle » 
(Fig. [T]). Dans le folklore de la chimie théorique, il s'agit des points de Pauling car Linus Pauling (prix 
Nobel : chimie, 1954; paix, 1962) avait le talent de savoir s'arrêter à ces points [3]. 




FIGURE 1 - Evolution du niveau de désaccord théorie-expérience en fonction de la sophistication du modèle théo- 
rique. Les points de meilleur accord sont les points de Pauling. Dans le cas des complexes à transition de spin, nous 
distinguons trois points de Pauling (voir texte). 

Dans le contexte des complexes à transition de spin, nous distinguons trois points de Pauling. Le premier 
est donné par la théorie du champ de ligands (LFT pour l'anglais Ligand Field Theory). Cette théorie 
semi-empirique et relativement ancienne est universellement utilisée pour comprendre le comportement 
des complexes de métaux de transition [3]. Il s'agit d'un véritable language, employé notamment par les 
expérimentateurs pour discuter leurs résultats en des termes tels que ceux de force du champ des ligands, de 
paramètres de répulsion électronique dits de Racah, ou encore d'effet néphélexautique. Cependant, la LFT 
a des limitations fortes. Ainsi, elle ne peut donner que très peu d'indications quantitatives sur la structure 
d'un complexe car son application nécessite de faire des suppositions sur la symétrie et donc sur la structure 
du complexe. Il est donc très difficile de prendre en compte de manière détaillée des effets dûs à de subtiles 
variations dans l'arrangement spatial des ligands autour du métal, ou d'inclure les effets d'environnement. 
Malgré ces difficultés, la LFT reste incontournable en tant que dépositaire de décennies d'expériences. 

Le troisième point de Pauling est donné par les calculs ab initio. Le terme latin ab initio veut dire « dès 
le début » et signifie dans notre contexte « à partir des principes premiers». Il désigne ici les calculs basés 
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le plus souvent sur un seul déterminant de Slater (Hartree-Fock) ou sur plusieurs déterminants de Slater et 
faits de manière auto cohérente, comme par exemple la méthode CASSCF (pour l'anglais complète active 
space self- consistent field) [5]. D'autres méthodes ab initio utilisés dans les études des complexes à transition 
de spin sont les méthodes CASPT2 [CASSCF avec correction perturbative au deuxième ordre) [6j [7] et 
SORCI (pour l'anglais spectroscopy oriented configuration interaction) jHlllO]. Leur avantage principal est 
une excellente capacité prédictive pour les petites molécules, des systèmes pour lesquelles elles servent de 
points de repère pour des méthodes plus simples comme la DFT. Leur désavantage est que les ressources 
informatiques requises (leur « coût de calcul») augmentent très rapidement avec la taille de la molécule. 
Cela a limité la fiabilité des études CASSCF pionnières des complexes à transition de spin [II]. Suite 
aux améliorations récentes dans la méthodologie CASSCF, des calculs ab initio bien plus fiables sur ces 
complexes commencent à apparaître |12H14| , Cependant, le pessismisme reste de mise pour les molécules 
d'au delà une certaine taille (~100 atomes). De plus, pour les petites molécules, une précision chimique 
supérieure a 5 kcalmol" 1 (1750 cm 1 ) reste très difficile à obtenir dès qu'un atome de métal de transition 
est présent. Nous appelerons ce niveau de précision, la meilleure précision anticipée pour les calculs sur 
les systèmes à métaux de transition (MPA-SMT). Pour bien interpréter ce qui suit, il est utile de garder 
à l'esprit ce chiffre de MPA-SMT = 1750 cm" 1 ! Pour atteindre une telle précision, certaines méthodes ab 
initio poussées telle que la CASPT2 ont recours elles aussi à un peu de semi-empirisme |15] , 

Le deuxième point de Pauling est la méthode qui probablement approche actuellement le plus du « juste 
milieu » d'Aristote. Il s'agit de la théorie de la fonctionnelle de la densité (DFT ; John Pople et Walter 
Kohn, prix Nobel de chimie de 1998). La DFT a des prédécesseurs importants |16H19j . Son formalisme 
moderne est basé sur deux théorèmes de Hohenberg et Kohn [20j dont le premier établit que l'Hamiltonien 
moléculaire est déterminé à une constante additive près par la densité de charge p de l'état fondamental. 
Le second établit que l'énergie et la densité de l'état fondamental peuvent être déterminées en minimisant 
une fonctionnelle (c'est-à-dire, une fonction d'une fonction) de la densité. Cette fonctionnelle est connue 
formellement |2T], mais est difficile à approximer en pratique. C'est la raison pour laquelle la formulation 
dominante de la DFT est celle de Kohn et Sham |22| dans laquelle la densité est exprimée à l'aide d'orbitales 
moléculaires (OMs) obtenues par résolution autocohérente des équations Kohn-Sham. Dans la formulation 
Kohn-Sham de la DFT, la quantité clé est l'énergie d'échange et de corrélation (énergie- :rc), E xc [p], qui 
est une fonctionnelle de la densité de charge p et regroupe les termes non-classiques intervenant dans la 
décomposition Kohn-Sham de l'énergie du système électronique. Mais, comme on ne connait pas de forme 
à la fois exacte et pratique de la fonctionnelle E xc , il faut l'approximer. 

Le point de départ de toute approximation de la fonctionnelle E xc est le gaz d'électrons homogène pour 
lequel des formes analytiques très précises de l'énergie- rcc par particule ont été établies. Si on suppose que la 
densité d'énergie-xc d'un système électronique réel est localement donnée par celle du gaz homogène d'élec- 
trons de densité celle du système au point considéré, on arrive à l'approximation locale dite LDA (pour 
l'anglais, local density approximation). Cette approximation simple marche étonnement bien pour l'optimi- 
sation des géométries moléculaires et le calcul des fréquences de vibration. Mais elle a plûtot tendance à 
surestimer les énergies de liaison et donc sous-estimer leurs longueurs. Pour aller au delà de la LDA, on 
ajoute une dépendance vis-à-vis du gradient de la densité pour arriver à l'aproximation GGA (de l'anglais, 
generalized gradient approximation) ; celle-ci donne des énergies de liaison bien meilleures que celles LDA. 
Malheureusement, il est difficile de faire mieux dans le cadre pure du formalisme Kohn-Sham dans lequel 
la densité électronique intervient comme seule « variable ». La tendance actuelle est donc d'étendre la dé- 
pendance de la fonctionnelle-xc aux orbitales Kohn-Sham, qui sont en fait des fonctionnelles implicites de 
la densité. Perdew a résumé la hiérarchie qui résulte des divers degrés de sophistication des fonctionnelles 
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TABLE 1 - L'échelle de Jacob des fonctionnelles [531 121] • (Une version plus à jour est donnée dans la Réf. |25|.) 



Le ciel de la chimie quantique 



double-hybride 




p(r),x(r),T(r),4> i (r),ip a (r) h 


hybride 




p(r),x(r),r(r),-0 î (r) 9 


mGGA c 




p(r),x(r),r(r) e , V 2 p(r)' 


GGA b 




p(r),x(r) d 


LDA a 







Le monde terrestre du modèle Hartree 



a LDA, voire texte. b GGA, voire texte. c Meta G G A. d Le gradient réduit x(r) = |Vp(r)|/p 4 / 3 (r). e L'énergie 
cinétique locale r(r) = ^2 p npipp(r)X7 2 îp p (r). f II y a des indications que la densité locale d'énergie cinétique, r(r), 
et le Laplacien de la densité de charge, V 2 /?(r), comprennent des informations semblables fJU] . 9 Orbitales occupées. 
h Orbitales vacantes. 

en une échelle de Jacob |23[ [25] montrée dans la Table [T] En montant cette échelle, on a plus de flexibilité 
pour approximer la fonctionnelle- xc et donc plus de possibilité d'améliorer l'accord expérience-théorie. Mais 
cette amélioration n'est pas assurée. Par contre, les calculs deviennent systématiquement plus coûteux. // 
est donc également important de pouvoir descendre l'échelle de Jacob. Pour plus d'informations sur la DFT, 
nous recommendons, par ordre de difficulté croissante, les livres |27H29j . 

La DFT est par sa nature limitée à l'état fondamental et à l'étude des propriétés indépendantes du 
temps. Nous pouvons tout au plus espérer qu'elle s'applique à l'état le plus bas en énergie de chaque symétrie 
d'espace et de spin, bien que les fonctionnelles usuelles n'aient pas la dépendence en symétrie que devrait 
avoir, tout au moins théoriquement, la fonctionnelle- xc |30H32j . La DFT dépendante du temps ou TDDFT 
est justement formulée pour traiter les problèmes dépendant du temps et les états excités. Elle repose sur 
des théorèmes semblables aux théorèmes de Hohenberg et Kohn et admet une formulation rappelant celle 
Kohn-Sham de la DFT traditionnelle |33j . De nouveau, une fonctionnelle- xc inconnue doit être approximée. 
Il est généralement admis que le champ autocohérent intervenant dans l'équation Kohn-Sham dépendant 
du temps réagit instantanément et sans mémoire à chaque variation temporelle de la densité. Il s'agit de 
l'approximation adiabatique qui permet d'utiliser pour les calculs TDDFT les fonctionnelles développées 
dans le cadre de la DFT. La caractérisation des états excités en TDDFT repose sur l'étude de la réponse de 
l'état fondamental à une perturbation dépendante du temps. Cette étude mène dans le cadre de la théorie 
des réponses au système d'équations, 

A(wj) B(wj) 

dont la résolution donne les spectres d'absorption électronique des molécules [34]. ojj est l'énergie d'excita- 
tion électronique verticale, et la force d'oscillateur associée, //, est calculée à partir des vecteurs Xj et Yj. 
Pour plus d'informations sur la TDDFT, on pourra se référer aux revues |24[ [34Tl38j ou aux livres |39[ |4"0] . 

Finalement, ajoutons qu'un des défis liés à la modélisation des systèmes à transition de spin est celui de 
traiter des systèmes de plus en plus grands. La méthode idéale serait à la fois simple et précise, mais avec 
une difficulté n'augmentant que lentement avec la taille du système. De ce point de vue, la DFT se prête 
a une reformulation dont la difficulté de calcul ne monte que proportionnellement avec le nombre d'atomes 
du système (méthodes d'ordre- N) [H] ; et il en est de même pour la TDDFT. 




1 

-1 




(1) 
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3 Application de la DFT à l'étude de la structure et de l'énergétique des 
complexes de métaux de transition dans les états bas-spin et haut-spin 



Nous discutons dans cette section la performance de la DFT pour la description des propriétés structurales 
et énergétiques des complexes de métaux de transition dans les états bas-spin (BS) et haut-spin (HS). Nous 
nous limitons pour cela aux complexes hexacoordinés d e de fer (II). L'état électronique fondamental d'un 
tel complexe de symétrie octaédrique (Oh) est, suivant la force du champ des ligands, l'état BS singulet 
1 Ai 9 (t2 9 ) caractérisé par un appariement maximum des électrons des niveaux d (ligands à champ fort), 
ou l'état HS quintuplet 5r ^2g(t2g e g) caractérisé par un appariement minimum des électrons des niveaux 
d et une multiplicité de spin maximale (ligands à champ faible). La dégénérescence orbitale de l'état HS 
est levée par des distorsions structurales dues à l'effet Jahn-Teller ; mais celles-ci sont en général faibles 
pour ces complexes d e octaédriques dans l'état HS. Nous allons d'abord examiner la performance de la 
DFT pour la détermination des géometries BS et HS des complexes, puis sa performance pour le calcul 
des différences d'énergie HS-BS. Nous verrons alors que, si la question de savoir si la performance des 
fonctionnelles s'améliore avec le temps est générale et importante, elle s'avère critique en ce qui concerne 
l'énergétique des états de spin des complexes. 

3.1 Performance de la DFT pour le calcul des géométries bas-spin et haut-spin 

La DFT donne en général une bonne description des géométries BS et HS des complexes de métaux de 
transition. Afin de discuter cela, nous prenons l'exemple du complexe [Fe(H20)g] 2+ - H s'agit d'un complexe 
dit HS car l'état HS est son état électronique fondamental et que celui-ci reste le seul peuplé à toutes 
températures. Les optimisations de géométries effectuées avec différentes fonctionnelles et sans contrainte 
de symétrie pour le complexe dans l'état HS mènent toutes à un état électronique quintuplet non dégénéré 
et à des géométries qui ne sont que légèrement distordues par l'effet Jahn-Teller. Pour [Fe(H20)g] 2+ dans 
l'état BS, ces optimisations mènent à un état singulet non dégénéré et à des géométries de symétrie Oh- Pour 
décrire ces structures, nous ne considérons que la longueur de liaison Fe-0 moyenne dont nous donnons les 
valeurs HS (rns) et BS (t'es) dans la Table[2] [Fe(H20)e] 2+ étant une espèce HS, seule sa structure HS a été 
caractérisée expérimentalement, en solution ou dans des composés cristallins. Les valeurs expérimentales 
de thSî données également dans la Table [21 varient de presque 0.09 Â : cette très grande dispersion souligne 
l'influence importante que peut avoir l'environnement sur les propriétés, ici structurales, des complexes. 
Cet allongement important de la liaison métal-ligand suite au changement d'états de spin, BS — > HS, est 
due à la promotion de deux électrons des orbitales tig non liantes aux orbitales e g anti-liantes. 

L'absence de données expérimentales pour [Fe(H 2 0) 6 ] 2+ dans l'état excité BS et le fait que les données 
expérimentales sont pour les molécules en phases condensées et que nos calculs sont effectués en phase 
gaz font qu'il est difficile de savoir quelles sont les valeurs prédites de tes, ?"hs et ArnB les plus fiables. 
Nous pouvons néanmoins affirmer que toutes les valeurs rns calculées sont en bon accord avec l'expérience, 
surtout si l'on prend en compte la probable compression de la structure de la molécule en phase condensée. 
Cet accord confirme que, comme anticipé, une force de la DFT est la qualité des géométries optimisées, ce 
qui permet d'aller bien au delà de la simple LFT. Nous pouvons également comparer les performances des 
différents types de fonctionnelles. La LDA tend à surestimer la force des liaisons. Cela se confirme ici par le 
fait que, comparée aux autres fonctionnelles, elle donne des liaisons Fe-0 trop courtes. Conformément à leur 
but historique, les GGAs corrigent cette surestimation de la LDA et donnent des liaisons plus longues. Les 
fonctionnelles hybrides incorporent une contribution d'échange HF à une forme GGA. Comme illustré ici 
avec la fonctionnelle B3LYP, elles montrent pour la détermination des structures des complexes de métaux 
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TABLE 2 - Comparaison théorie-expérience des distances (Â) moyennes Fe-0 pour le complexe [Fe(H20)e] 2+ - La 
notation utilisée est « méthode de calcul/base utilisée » (le nom d'une base est le sigle utilisé dans l'article cité). 
La notation utilisée pour les fonctionnelles est : x (échange pour l'anglais exchange), c (corrélation), xc (échange- 
corrélation). Souvent, l'abbréviation d'une fonctionnelle est composée du sigle de la partie d'échange suivi du sigle de 
la partie de corrélation (e.g., BP86 est composée de la GGA-x B plus la GGA-c P86). Pour les autres abbréviations, 
se reporter à la Table Q] et à la discussion associée en Sec. [5] 



Méthode 






r HS 




LDA/A m :1965 A (LDA-x) [H] 




1.912 


2.070 


0.158 


BP86/E a : 1988 (GGA-x) g2] + 


1986 (GGA-c) g3] 


1.996 


2.145 


0.149 


BLYP/E Q :1988 (GGA-x) g2] + 


1988 (GGA-c) gj 


2.025 


2.166 


0.141 


PW91/E a : 1991 (GGA-xc) g5J B 


1.992 


2.137 


0.145 


PBE/C" a : 1996 (GGA-xc) g5] 




2.011 


2.154 


0.143 


RPBE/C" Q : 1999 (GGA-x) g7] 


+ 1996 (GGA-c) [1 


3 2.048 


2.187 


0.139 


B3LYP/E a : 1994 (Hybride-xc) c 




2.018 


2.154 


0.136 


Hartree-Fock b 




2.080 


2.193 


0.113 


CASSCF(10,12)/II b 




1.989 


2.113 


0.124 


CASSCF(12,10) a 




2.077 


2.186 


0.119 


Exp. a 






2.068-2.156 




LFT 








0.126 



A Paramétrisation VWN de la LDA-c gS]. B La fonctionnelle PW91 était annoncée en congrès et était rapidement 
propagée dans le monde DFT sous la forme d'un morceau de code FORTRAN, mais pendant longtemps il manquait 
une publication officielle décrivant la fonctionnelle. c La fonctionnelle B3LYP est identique à une fonctionnelle 
hybride proposée par Becke gH], sauf que la composante GGA-c, P91 gS], est arbitrairement remplacée par la 
GGA-c LYP gj. a Résultats tabulés dans la Réf. [50]. b Réf. [H]. c L'estimation LFT est décrite en Réf. [5T]. 

de transition des performances similaires à celles des GGAs. 

Bien que le choix de la fonctionnelle n'ait pas d'influence dramatique sur les géométries HS et BS 
calculées, il affecte de manière remarquable la différence ArnB) dont les variations approchent les 40%. 
De façon générale, il reste difficile de déterminer quelle fonctionnelle est la meilleure pour le calcul des 
géométries, soit parce que les effets d'environnement ne sont pas pris en compte dans les calculs, soit parce 
que la taille du système empêche d'obtenir pour le système isolé une géométrie de référence à partir de 
méthodes ab initio poussées. Dans le cas du [Fe(H20)e] 2+ , en nous fiant à l'estimation LFT (basée sur 
des décennies d'expérience) de ArjjB et les calculs CASSCF, nous pouvons éliminer la méthode HF (où 
le traitement de la corrélation électronique est absente) et celle de la LDA au profit des fonctionnelles de 
type GGA ou hybride. Cette conclusion est valable pour la plupart des complexes à métaux de transition. 

3.2 Performance de la DFT pour l'énergétique des états de spin 

Par énergétique des états de spin, il nous faut distinguer la différence des énergies de point zéro HS et BS, 
AE^ B = E^ s - E BS , de ses contributions électronique, AE'hb = -^HS - ^BS 3 et vibrationnelle, AE^ B = E^ s 
- E BS . AE'hb est accessible par optimisation des géométries HS et BS, et AE^ B par calcul des fréquences 
du complexe dans les état BS et HS. Les calculs DFT donnent au niveau GGA ou hybride de très bonnes 
estimations des fréquences des complexes et, par conséquent, de AE^ B . La situation se complique lorsqu'il 
s'agit de déterminer A.Ehb qui, dans les complexes à transition de spin, prend des valeurs de l'ordre de 
100~1000 cm -1 . Pour de tels calculs, les meilleures méthodes ab initio ont une meilleure précision anticipée 
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MPA-SMT d'environ 1750 cm 1 (Sec. [2]). Le défi est donc bien grand pour la DFT '. 

La situation actuelle peut se résumer de la manière suivante. La fonctionnelle LDA tend à fortement 
surestimer A£hb- Les fonctionnelles GGAs corrigent les défauts de la LDA. Mais celles des études DFT 
d'avant le début des années 2000 surestiment toujours la stabilité de l'état BS par rapport à l'état HS. 
Ces mêmes études montrent qu'à l'inverse, l'hybride B3LYP, alors la fonctionnelle la plus sophistiquée 
couramment utilisée, surestime la stabilité de l'état HS par rapport à l'état BS : l'écart A£"hb est sous- 
estimé, comme avec la méthode HF. Le début des années 2000 a vu la performance des méthodes DFT 
s'améliorer. C'est le cas pour les hybrides avec la B3LYP* obtenue empiriquement en ajustant la contribution 
d'échange HF intervenant dans la B3LYP, de 20% à 15%, de façon à reproduire les états fondamentaux 
d'une série de complexes |52] , II est important de se permettre de descendre l'échelle de Jacob (Sec. [2]), et 
l'apparition de nouvelles GGAs plus sophistiquées en a donné l'occasion. Toutefois, lorsqu'il s'agit de tester 
les fonctionnelles pour leur capacité à prédire AI?hb> la plus grande difficulté se trouve être la collection et 
la sélection des données de comparaison. Cela est loin d'être trivial car les calculs de structures électroniques 
sont généralement effectués sur le complexe isolé, alors que les données expérimentales sont souvent pour 
des complexes en milieu condensé et dépendent plus étroitement de l'enthalpie ou enthalpie libre que de 
l'énergie électronique. A titre d'exemple, nous présentons dans le Tableau [3] des résultats obtenus pour les 
valeurs de AE R b dans les complexes [Fe(H 2 0) 6 ] 2+ , [Fe(NH 3 ) 6 ] 2+ et [Fe(bpy) 3 ] 2+ (bpy = 2,2'-bipyridine). 
Des repères de qualité ont été obtenus pour les petits complexes HS [Fe(H20)g] 2+ [50j et [Fe(NH3)g] 2+ 
par des calculs ab initio poussés CASPT2 et SORCI [501 EU- Quant au complexe BS [Fe(bpy)3] 2+ , une 
estimation fiable de A£"hb a été obtenue par l'étude expérimentale de la cinétique de relaxation HS —¥ BS 
du complexe dopé dans différentes matrices |53| 15^] . 

L'incertitude sur les résultats des calculs ab initio est 1013 cm -1 pour le complexe [Fe(H20)e] 2+ et 
2105 cm" 1 pour le complexe [Fe(NH3)g] 2+ , ce qui est cohérent avec la MPA-SMT de 1750 cm -1 . L'incerti- 
tude sur la valeur expérimentale trouvée pour le complexe [Fe(bpy)3] 2+ est plus grande (2500 cm -1 ) ; cela 
reflète l'influence du milieu sur AE^hb- Comme cela pouvait être anticipé, la LDA et les GGAs d'avant 1998 
(BP86, PW91, RPBE) surestiment la différence A£"hb- On note également que, pour [Fe(bpy)3] 2+ , AEhb 
est sous-estimée par les hybrides classiques B3LYP et PBE0. La situation est cependant moins simple dans 
le cas du complexe [Fe(H20)g] 2+ pour lequel A£"hb est sous-estimée par la PBE0 mais est raisonnablement 
approchée par la B3LYP. Enfin, les deux fonctionnelles hybrides surestiment toutes deux la meilleure esti- 
mation de AEhb dans le cas du [Fe(NH3)g] 2+ . La fonctionnelle reparamétrée B3LYP* amène par rapport à 
la B3LYP une nette amélioration pour l'évaluation de A£"hb dans le complexe [Fe(bpy)3] 2+ . Par contre, elle 
ne fait qu'accentuer la tendance de la B3LYP à surestimer AE'hb dans le cas des complexes [Fe(H20)e] 2+ 
et [Fe(NH3)g] 2+ qui sont chimiquement très différents de ceux utilisés pour sa paramétrisation. Les GGAs 
post-1996 telles que la RPBE, les "HCTH" et les fonctionnelles OLYP et OPBE donnent des résultats 
bien meilleurs que les GGAs précédentes. Les GGAs OLYP et OPBE ne diffèrent que par la fonctionnelle 
de corrélation dont le choix s'avère moins important pour l'énergétique des états de spin que celui de la 
fonctionnelle d'échange. Ces deux GGAs sont actuellement deux des fonctionnelles les plus performantes 
pour l'énergétique des états de spin des complexes de fer(II). Ainsi, la OLYP donne de meilleurs résultats 
que la B3LYP* pour le calcul des différences d'enthalpie libre AGhb de la série de complexes utilisée pour 
paramétrer la B3LYP* |59] , Il semblerait qu'en fait la contribution d'échange HF dans la B3LYP* doive être 
abaissée en dessous de sa valeur de 15% pour améliorer la performance de cette hybride |59[ |6T] , La dépen- 
dance de A.Ehb vis-à-vis de la contribution d'échange HF a aussi été étudiée dans le cas du [Fe(bpy)3] 2+ ; 
les résultats suggèrent qu'effectivement cette contribution doit être réduite à environ 10% aussi bien dans 
la B3LYP que dans l'hybride PBE0 |53] , Ajoutons qu'une étude récente a montré la bonne performance 
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TABLE 3 - Comparaison des fonctionnelles pour les différences énergétiques A£"hb = -Ehs — -Ebs (cm -1 ) dans 
quelques complexes de fer(II) bien caractérisés. Le spin de l'état fondamental est entre parenthèses. Pour le com- 
plexe Z?3 [Fe(bpy)3] 2+ (bpy = 2,2'-bipyridine), nous supposons que l'état HS est la composante trigonale 5 Ai de 
l'état champ-des-ligands 5 T2g ; toutefois les deux composantes trigonales 5 E et 5 Ai sont quasi-dégénérées [53]. Les 
fonctionnelles sont en ordre chronologique. 



Méthode Molécule 
fonctionnelle : année (type) [référence] [Fe(H 2 0) 6 ] 2+ (HS) [Fe(NH 3 ) 6 ] 2+ (HS) [Fe(bpy) 3 ] 2+ (BS) 

LDA : 1965 A (Locale- s) [22] -3 896° 8 187 b 
BP86 : 1988 (GGA-x) 

+ 1986 {GGA-c) [43] -8 798 a 790 b 11 135 d 

PW91 : 1991 (GGA-xc) [45] B -9 232° 617 b 11 887 d 

B3LYP : 1994 (Hybride-zc) [— c \ -11 465 a -4 978 b 1 211 d 

PBE : 1996 (GGA-xc) 05] -10 081° -147 b 11 337 d 

HCTH93 : 1998 (GGA-xc) [55] -18 779 b -9 430 b 

HCTH147 : 1998 (GGA-xc) [55] -18 211 b -8 576 b 

HCTH407 : 1998 (GGA-xc) [55] -19 631 b -10 682 b 

VSXC : 1998 (Méta- GGA-xc) [5J] -13 975 e ' / -5 928^ -1 822 e ^ 

RPBE : 1999 (GGA-xc) 07] -11 844 a -2 744 b 6 640 d 

PBEO 13 : 1999 (Hybride- xc) [57] -14 504 b -7 195 b -466 d 
OLYP B : 2001 (GGA-x) [55] 

+ 1988 (GGA-c) [S] -16 284 e ^ -6 2 684 e ^ 

OPBE B : 2001 (GGA-x) [55] 

+ 1988 (GGA-c) gB] -17 200 9 -6 650 9 5 210 9 

B3LYP* : 2002 (Hybride- xc) [5J -10 456 e < / -3 226^ 3 332 e ^ 

Meilleures estimations 

Exp. 3 500 — 6 000 d 

Ab initio c -13 360 — -12 347 a -11 230 — -9 125 b 

Moyenne de l'intervale -12 853 -10 177 4 750 



A c Se référer à la Table [U D Aussi connue par le nom PBE1PBE. E La GGA-x OPTX [55] est souvent désignée 
uniquement par la lettre O, comme dans les fonctionnelle-xc OLYP et OPBE. a Réf. [50]. b Réf. [5J]. c Résultats 
de calculs CASPT2 et SORCI avec, pour les résultats CASPT2, une correction atomique empirique [501 EU- d 
Réf. [53]. e Réf. [5J5]. ' A la géométrie optimisée avec la même fonctionnelle, mais pour ces complexes, le choix de la 
fonctionnelle influe peu sur les résultats des optimisations de géométries. 9 Réf. [60J. 



de la OPBE pour l'énergétique des états de spin de complexes de fer (II) et aussi de fer(III) |60j . 

// apparaît donc qu'en passant de l'hybride B3LYP* aux GGAs OLYP et OPBE, nous avons réussi à 
descendre l'échelle de Jacob sans perte notable de précision. Cela permet d'envisager avec confiance l'étude 
de systèmes de métaux de transition beaucoup plus grands que ceux étudiés jusqu'ici. Nous pouvons aussi 
raisonnablement espérer une amélioration régulière de la performance de la DFT pour l'énergétique des états 
de spin. Le domaine ne peut en effet que bénificier du développement constant de nouvelles fonctionnelles, de 
degrés de sophistication plus ou moins élevés, et améliorant la description DFT des propriétés moléculaires. 

Une difficulté majeure liée à l'utilisation de la DFT pour l'énergétique des états de spin est qu'il 
est difficile et risqué de préjuger de la performance d'une fonctionnelle pour des complexes sur la base 
de sa bonne, voire excellente performance pour d'autres complexes. Il serait donc intéressant de pouvoir 
estimer Ai?HB en contournant cette dépendance vis-à-vis de la fonctionnelle. Ainsi, en admettant que 
les différentes fonctionnelles- xc donnent des surfaces d'énergie potentielle de formes très semblables mais 
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décalées en énergie les unes des autres, il a été montré [62] que la différence A(A.Ehb) des différences 

d'énergie A£"hb calculées pour deux complexes Ci et C2, A(A£hb) = A.Ehb(Ci) — A^hb(C2), ne devrait 
que peu dépendre de la fonctionnelle utilisée. Cette proposition a été vérifiée à partir de calculs effectués sur 
deux complexes de fer(II) à transition de spin en utilisant 53 fonctionnelles dont des GGAs, des hybrides 
et des méta-GGAs |62] . Les GGAs utilisées ont donné de très bons résultats, les hybrides ont fait un 
peu moins bien, et certaines méta-GCMs échouent tout simplement. Néanmoins, il apparaît donc possible, 
connaissant précisément AE'hb pour un complexe de référence, d'estimer par calculs DFT sa valeur dans 
d'autres complexes avec une précision semblable, ou éventuellement supérieure, à la MPA-SMT. 



Parmi les méthodes capables de décrire les états excités des molécules, la TDDFT est très probablement celle 
la plus utilisée actuellement pour l'étude des molécules de taille moyenne et grande. Par souci de concision, 
nous nous limitons ici à quelques observations sur la place de la TDDFT dans l'étude des complexes de 
métaux de transition, et plus particulièrement ceux à transition de spin. 

Un domaine d'applications très important de la TDDFT est la modélisation des propriétés photophy- 
siques des complexes de ruthénium(II), des complexes cl® à couche fermée dont l'intérêt vient notamment 
de leur utilisation comme pigments photosensibles dans les cellules photovoltaïques de Grâtzel [63J. Dans ce 
domaine, le nombre de travaux basés sur les calculs TDDFT a véritablement explosé du fait de la disponi- 
bilité de la méthode dans la majorité des logiciels de chimie ou physique quantique, ainsi que du bon accord 
expérience-théorie auquel on peut aboutir. En fait, pour ce qui est de l'application de la TDDFT à l'étude 
des molécules à couches fermées comme ces complexes de ruthénium(II), la situation est relativement simple 
car l'on doit essentiellement se préoccuper du choix de la fonctionnelle à utiliser. De plus, la plupart des 
fonctionnelles GGAs et hybrides permettent une bonne voire très bonne description des états excités de 
valence de ces molécules. Cependant, du fait de la nature approximée des fonctionnelles et de l'approxima- 
tion adiabatique qui sous-tend son application usuelle (Sec. [4]), la TDDFT connaît des limitations. Ainsi, la 
caractérisation des états excités hauts en énergie (états de Rydberg) requiert l'utilisation de fonctionnelles 
avec un comportement asymptotique correct. De même, des difficultés existent lorsqu'il s'agit de décrire les 
états excités ayant un caractère de transfert de charge à longue distance ou mettant en jeu des excitations 
électroniques multiples. Nous renvoyons le lecteur aux Refs. |24[ [33HJQ] pour des discussions détaillées des 
limitations actuelles de la TDDFT et des développements auxquels elles donnent lieu. 

La TDDFT est également appliquée à l'étude des propriétés optiques d'espèces à couche ouverte telles 
que le complexe cf [Co(bpy)3] 2+ . Les spectres d'absorption et de dichroïsme circulaire de ce complexe chiral 
ont ainsi été modélisés par des calculs TDB3LYP/TZVP |64| . La Fig. [2] montre les spectres obtenus pour 
l'énantiomère A du complexe dans l'état BS. Les forces d'oscillateur et les forces rotatoires (voir ci-dessous) 
calculées sont également représentées par des diagrammes en bâtons. Pour simuler le spectre d'absorption 
e = e{y) (e en M -1 cm -1 , v en cm -1 ), le profil d'absorption ei{p) associé à la I-ème transition centrée en vj 
a été calculée en convoluant la force d'oscillateur // par une Gaussienne de largeur à mi-hauteur 2000 cm -1 
dont la normalisation garantit l'égalité 



Le pouvoir rotatoire Rj de cette transition est liée à la différence Ae/ entre les coefficients d'absorption 



4 États excités par TDDFT 




(2) 
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FIGURE 2 - Spectres d'absorption (gauche) et de dichroïsme circulaire (droite) calculés pour l'énantiomère A 
du complexe [Co(bpy)3] 2+ dans l'état BS (calculs TDB3LYP/TZVP). Les spectres d'absorption et de dichroïsme 
circulaire mesurés sur un monocristal du composé [Co(bpy)3][LiRh(ox)3] à 11 K et 15 K, respectivement, sont 
également représentés. D'après la Réf. |64| . 

pour des lumières circulaires polarisées gauche (e s ) et droite (e^) (Ae = e g — e r ), par la relation 

Ri = (22.97 x 10~ 40 cgs M cm) / — ^ dû . (3) 

Jo v 

Pour simuler le spectre de dichroïsme circulaire, Ae/(P) a été calculé en convoluant Rj par une Gaussienne 
de largeur à mi-hauteur 2000 cm -1 dont la normalisation garantit l'égalité ci-dessus. La Fig. [2] permet 
la comparaison des spectres calculés avec ceux mesurés sur un monocristal du composé à transition de 
spin [Co(bpy)3][LiRh(ox)3] qui, aux températures cryogéniques des mesures, est dans l'état BS. L'accord 
théorie-expérience est très bon : [Co(bpy)3] 2+ a donc la configuration A dans le monocristal étudié [64] . 

Bien que la TDDFT permette généralement de prédire l'enveloppe correcte d'un spectre, elle ne permet 
pas nécessairement de bien résoudre celui-ci en donnant une bonne description des états impliqués. Cela 
dépend de la qualité de la fonctionnelle et, pour les systèmes à couche ouverte, du problème de contamination 
de spin dans l'état fondamental et dans les états excités. La contamination de spin se produit lorsque la 
fonction d'onde approximative ^ d'un état de spin S se trouve ne pas être un état propre de l'opérateur 
de spin S 2 de valeur propre S(S + 1). Elle est diagnostiquée par l'écart de la valeur moyenne calculée 
(S 2 ) = |S 2 | *) à celle attendue S (S + 1). Ce problème est moins important dans les calculs DFT que 
dans ceux Hartree-Fock. Dans les calculs TDDFT, son absence complète (ou partielle) est une condition 
nécessaire, bien que non suffisante, pour une description correcte de l'état excité. Grâce aux récents travaux 
des Réfs. |24 [ I65 [ [66] . nous avons désormais des outils pour assigner automatiquent le spin des états excités 
en TDDFT, et ainsi identifier et traiter une éventuelle contamination de spin. La Fig. [3] illustre cette 
approche dans le cas du complexe [Fe(H20)e] 2+ [67] • Elle présente les courbes d'énergie potentielle des 
états quintuplets fondamental (Qo) e t excités (Qj>i) du complexe [Fe(H20)e] 2+ le long de la coordonnée 
de distance Fe-0 associée au mode de respiration totalement symétrique du complexe. Les évolutions des 
valeurs moyennes de (S 2 ) calculées pour ces états le long des courbes d'énergie potentielle y sont également 
représentées. Les quatres courbes les plus basses en énergie correspondent à l'état fondamental et à des 
états excités de type d — » d. Ces états sont libres de contamination de spin ((S 2 ) = 6). Les courbes plus 
hautes en énergie correspondent à des excitations de type transfert de charge métal — > ligand. Selon la 
géométrie du complexe, ces états présentent une contamination de spin ((S 2 ) = 7) et doivent alors être 
écartés ou corrigés pour leur rendre un sens physique [66j. 

La possibilité de contamination de spin est en partie une conséquence de la simplicité de la TDDFT dans 
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FIGURE 3 - Courbes d'énergie potentielle des états quintuplets (Qo) et excités (Qi>\) du complexe [Fe(H20)e] 2+ le 
long de la coordonnée de distance Fe-0 (bas), et valeurs moyennes de (S 2 ) calculées pour ces états le long des courbes 
d'énergie potentielle (haut). Calculs DFT au niveau PBE pour l'état fondamental et calculs TDDFT&u niveau LDA 
TDA pout les états excités. 



son état de développement actuel. Mais cette même simplicité permet une application très intéressante : 
la photodynamique semiclassique avec saut de surfaces qui est une technique de simulations de dynamique 
non-adiabatique |68H73] . L'étude de la dynamique nucléaire est souvent requise pour bien comprendre 
les rendements des réactions photo chimiques qui, pour la plupart, passent par des intersections coniques 
ou à leur proximité. De ce point de vue, l'étude de la photochimie de l'oxirane a été très instructive. Une 
première étude TDDFT de l'oxirane a montré l'importance en TDDFT conventionnelle^ de l'approximation 
Tamm-Dancoff [TDA pour l'anglais Tamm-Dancoff approximation) [B=0 dans l'Eq. ([T|] pour réduire les 
effets nuisibles dûs aux instabilités de spin dans les calculs TDDFT des surfaces d'énergie potentielle |68| . 
L'application de la photodynamique TDDFT semiclassique dans la TDA a ensuite donné des résultats 
en très bon accord avec le mécanisme classique de Gomer et Noyés qu'ils ont permis de compléter |71] , 
Ces bons résultats peuvent toutefois surprendre car il ne peut pas exister dans la TDDFT conventionnelle 
de vraies intersections coniques entre l'état fondamental et les états excités [74] . Les surfaces d'énergie 
potentielle des états Sq et Si obtenues par des calculs TDDFT TDA et CASSCF ont en fait la même allure 
au voisinage de l'intersection conique déterminée par les calculs CASSCF. La différence est qu'en TDDFT 
TDA, au lieu d'avoir une vraie intersection conique, les deux cônes s'interpénétrent. Ce détail est important 
pour déterminer la vitesse de saut de surface. Mais, comme les résultats de la Réf. |7T] le démontrent, il ne 
l'est pas forcément pour prédire le mécanisme de la photoréaction. 



x En TDDFT conventionnelle, la dépendance en fréquence des matrices A et B de l'Éq. JTJ est négligée. 
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5 Incorporation des effets environnementaux 



Nous nous sommes jusqu'ici uniquement intéressés aux propriétés des complexes isolés. L'influence de leur 
environnement est cependant loin d'être anodine. Prenons les complexes hexacoordinés de fer(II). Lors 
d'un changement d'états de spin BS — > HS, leur volume moléculaire augmente considérablement suite à 
l'allongement des longeurs des liaisons métal-ligand (Table [2]). Cette grande variation de volume est à 
l'origine de la sensibilité du phénomène de transition de spin à la pression externe et à des modifications 
dans la seconde sphère de coordination des complexes |54| . En effet, l'application d'une pression entraîne 
une déstabilisation de l'état HS par rapport à l'état BS de plus faible volume moléculaire, et donc une 
augmentation de la différence des énergies de point zéro, A££ B . Pour le complexe BS [Fe(bpy) 3 ] 2+ , AE^ B 
prend des valeurs comprises entre 2500 et 5000 cm" 1 suivant l'environnement du complexe |53[ 175] , Par 
analogie, l'influence de l'environnement sur AE^ B peut aussi être rationalisée en termes d'une « pression 
interne » ou « pression chimique » qui déstabilise l'état HS par rapport à l'état BS. Toutefois, cela ne 
permet pas d'accéder aux détails des interactions hôte-invité mises en jeu. Une analyse approfondie de ces 
interactions est cependant possible par calculs DFT, comme cela a été récemment montré pour le composé 
d'inclusion [Fe(bpy)3] 2+ @Y à partir du modèle supramoléculaire de la Fig. |4] [76J. 



FIGURE 4 - (Gauche) Modèle supramoléculaire de symétrie C3 du composé d'inclusion [Fe(bpy) 3 ] 2+ @Y : le modèle 
comprend 229 atomes et est constitué du complexe et des atomes environnants d'oxygène et de silicium qui définissent 
la supercage ; la valence des atomes de silicium est saturée par des atomes d'hydrogène qui ont été omis par souci de 
clarté. (Droite) Comparaison des géométries BS de [Fe(bpy)3] 2+ en phase gaz (noir) et incorporé dans la supercage 
de la zéolite Y (blanc) : l'encagement entraîne une légère contraction et distortion de la géométrie du complexe. 
Figure reproduite avec permission de la Réf. |76) . 

Les structures de [Fe(bpy)3] 2+ et du modèle de [Fe(bpy)3] 2+ @Y ont été optimisées dans les états BS et 
HS en utilisant différentes fonctionnelles G G As et hybrides. Pour chacun des systèmes, les structures BS ou 
HS obtenues avec les différentes fonctionnelles sont très semblables. Leur analyse montre que l'encagement 
entraîne une légère contraction et une légère distortion de la structure de [Fe(bpy)3] 2+ (Fig. HJ. L'influence 
de l'encagement sur la structure s'avère plus marquée dans l'état HS que dans l'état BS à cause du volume 
moléculaire plus grand du complexe dans l'état HS. L'éclatement quadrupolaire AEq du 57 Fe fournit une 
mesure de la distorsion du complexe et peut être déterminé expérimentalement par spectroscopie Môssbauer. 
Les valeurs de AEq obtenues par calculs DFT pour [Fe(bpy) 3 ] 2+ isolé et encagé ont pu être corrélées avec 
le degré de distortion prédit pour le complexe dans les états BS et HS (voir Réf. |76]). De plus, ces valeurs 
calculées sont en bon voire excellent accord avec les références expérimentales disponibles. 

La DFT permet d'arriver efficacement à une description très satisfaisante des structures et propriétés 
Môssbauer de [Fe(bpy)3] 2+ @Y dans les états HS et BS. Cependant, comme on pouvait l'anticiper (Sec. [3j, 
et comme le montre les données du Tab. [4j l'évaluation en DFT de l'énergétique des états de spin du 




12 



TABLE 4 - Valeurs calculées des différences d'énergie HS-BS AShb (cm -1 ) dans [Fe(bpy) 3 ] 2+ et [Fe(bpy) 3 ] 2+ @Y ; 
variation A(A-Ehb) de A£?hb suite à l'encagement ; et décomposition de A(A£'hb) en ses contributions géométrique 
(Ai^'g) et d'énergie d'interaction hôte-invité (A-Egg). D'après la Réf. [75]. 



PBE B3LYP* HCTH 03LYP OLYP 
Différence d'énergie HS-BS, AShb 

[Fe(bpy) 3 ] 2 + +10087 +3849 +141 -811 +3660 

[Fe(bpy) 3 ] 2 +@Y +12004 +5925 +2687 +1941 +6594 

A(A£ HB ) = AS HB [Y] - A£ HB [0] 

A(A£"hb) +1917 +2076 +2546 +2752 +2934 

A(A£ H b) = AE^ + A^b 

AE'^ +1246 +1267 +1621 +1286 +1331 

AE$g +671 +809 +925 +1466 +1593 



composé [Fe(bpy)3] 2+ @Y reste un problème difficile. En effet, la différence d'énergie AEhb[Y] calculée varie 
entre +1941 et +12004 cm -1 suivant la fonctionnelle. Néanmoins, toutes les fonctionnelles prédisent que 
l'encagement induit la déstabilisation de l'état HS par rapport à l'état BS. De plus, les valeurs trouvées pour 
la variation A(AE'hb) de A£"hb sont très cohérentes entre elles car comprises entre +1917 et +2934 cm -1 , 
et donc avec une variation inférieure à la MPA-SMT. La meilleure estimation de A(A£'hb) est donc [76] : 



A(A£ H b) = A£ H b[Y] -AE hb [0] 

= +2500 ± 1000 cm" 1 . (4) 

Les notations X[0] et AT [Y] désignent les valeurs de la grandeur X pour [Fe(bpy) 3 ] 2+ isolé et encaee, 
respectivement. Pour analyser l'influence des interactions hôte-invité sur l'énergétique des états de spin, 
A(AE'hb) est décomposée en deux contributions AE^'g et AE^. La première est d'origine structurale et 
s'écrit Ai?g|' = E^ — E^, où E^ lst est l'énergie requise pour déformer dans l'état de spin T le complexe 
et la supercage pris séparément de leurs structures relaxées à celles qu'ils ont dans [Fe(bpy)3] 2+ @Y. La 
seconde est la variation de l'énergie des interactions hôte-invité (E mt ) induite par le changement d'états de 
spin BS ^ HS : AE$ B = £gs ~ ^bs- 

Le Tab. |4]montre que la valeur de AE^ dépend peu de la fonctionnelle utilisée (AE^g ~ +1300 cm" 1 ), 
contrairement à celle de AE^g. La variation des valeurs de AE$g est donc la responsable principale de 
la dispersion des valeurs de A(AE'hb)- Elle est due à l'influence de la fonctionnelle sur la description des 
structures du complexe et de la supercage |76j . AE^ > car la distortion liée à l'encagement est plus 
grande dans l'état HS que dans l'état BS. On note également que AE^ > : les interactions hôte-invité 
sont donc moins stabilisantes suite au changement d'états de spin BS — > HS. L'énergie E mt des interactions 
instantanées entre deux fragments peut être décomposée suivant |77[ [75] 

^int ^elstat _|_ ^Pauli _|_ JTJ ^ (tj\ 

E cktat est l'énergie d'interaction électrostatique classique entre les distributions de charges non perturbées 
de l'hôte et de l'invité. E p 1 correspond à la répulsion de Pauli qui est responsable de toute répulsion 
stérique entre les deux entités. E ovh est l'énergie d'interaction orbitalaire, c'est-à-dire l'énergie gagnée en 
permettant à la densité électronique de se relaxer. Le Tab. [5]présentent les résultats obtenus par l'application 
de cette analyse à l'étude des interactions hôte-invité dans [Fe(bpy)3] 2+ @Y. Les interactions sont liantes 
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dans les deux états de spin (E^ < 0) et la liaison est plus électrostatique que covalente (pistât < £° rb < ; 
r = LS, BS). La contribution majoritaire à AE^ vient de la répulsion de Pauli, ce qui intuitivement 
s'explique par l'augmentation du volume du complexe. Mais AE^^ h ne représente que ~53% de AE^ : 
les interactions électrostatiques et orbitalaires sont moins stabilisantes et contribuent de manière comparable 
aux «47% restant de AE&L Cela est moins accessible à l'intuition. 



TABLE 5 - Analyse des interactions hôte-invité dans [Fe(bpy)3] 2+ @Y (énergies en cm 1 ; résultats OLYP) |76| . 







^elstat 


^Pauli 






r 


= BS 


-15558 


+ 17442 


-12158 


-10274 


r 


= HS 


-15251 


+18099 


-11874 


-9026 






y\ gjelstat 
HB 




A££ r B b 








+307 


+657 


+284 


+1248 



En fait, pour une géométrie donnée du système [M(bpy)3] 2+ @Y où M 2+ est un dication de métal de 
transition, les énergies d'interactions hôte-invité s'avèrent ne pas dépendre de l'état de spin du dication, ni 
même de sa nature. Cela a permis de définir les interactions hôte-invité dans le système [M(bpy)3] 2+ @Y 
comme étant les interactions à couches fermées entre la première sphère de coordination du complexe définie 
par les ligands et la deuxième sphère de coordination définie par la supercage, sous les effets polarisants du 
dication [76J. La bonne performance de la plupart des fonctionnelles modernes pour la description de ces 
interactions à couches fermées explique la bonne qualité des résultats obtenus. La meilleure estimation 
actuelle de AE UB [0] obtenue par calculs CASPT2 est de 3700 ± 1000 cm" 1 PJJ [13]. On obtient alors 
[Eq. (|4])] comme meilleure estimation pour l'énergétique des états de spin dans [Fe(bpy)3] 2+ @Y 

A£ H b[Y] = 6200 ± 1500 cm" 1 . (6) 

et une différence des énergies de zéro-point de A£^ B [Y] ~ 5400 cm" 1 [76J. Dans la théorie de la relaxation 
HS — > BS [75], cette valeur de AE'ggjY] correspond à une durée de vie ths = I/^hb^ 0) de l'état HS 
de moins de 10 ns. Cela est en accord avec la limite ths[Y] < 60 ns déduite d'expériences de spectroscopie 
Môssbauer d'émission |79] , 

L'étude du composé [Fe(bpy)3] 2+ @Y montre que la DFT permet de prendre en compte les effets d'en- 
vironnement de manière efficace et fiable et d'arriver à une analyse fine des interactions hôte-invité. Pour 
les composés cristallins, les calculs DFT périodiques peuvent être appliquées à l'étude approfondie des 
effets d'empilement et de contre-ions en phase solide. De façon similaire, l'étude des complexes en solution 
par simulations de dynamique moléculaire permettra de progresser dans la compréhension de l'influence 
du solvant et des fluctuations thermiques sur leurs propriétés structurales et électroniques. Ainsi, l'étude 
récente d'une solution aqueuse de [Fe(bpy)3](Cl)2 à 300 K par simulations de dynamique moléculaire |80] 
a permis de décrire les structures du complexe [Fe(bpy)3] 2+ et de sa couche de solvatation ainsi que leur 
dépendance vis-à-vis de l'état de spin (Fig. |5|). 

6 L'Avenir 

Des progrès majeurs ont été réalisés durant la dernière décennie dans la modélisation des complexes à 
transition de spin. En reprenant l'idée des trois points de Pauling présentée au début du chapitre, il nous 
faut d'abord préciser que la LFT (le premier point) n'a pas été à proprement parler supplantée ; celle-ci 
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FIGURE 5 - Etude par simulations de dynamique moléculaire Car-Parrinello d'une solution aqueuse de 
[Fe(bpy)3](Cl)2 à 300 K dans les états BS et HS : (gauche) cliché pris sur la trajectoire HS montrant [Fe(bpy)3] 2+ , 
un anion Cl - , des molécules d'eau et la densité de spin centrée sur le fer; (droite) cliché pris sur la trajectoire BS 
montrant la première sphère de solvatation du complexe qui est constituée de molécules d'eau intercalées dans les 
sillons définis par les espaces entre les ligands. D'après la Réf. |80| . 

restant particulièrement bien ancrée dans la communauté des expérimentateurs. Les avancées viennent du 
fait que, dans un effort concerté, les travaux de différents groupes ont convergé pour fermement établir 
le second point de Pauling (la DFT et la TDDFT). La DFT est non seulement performante pour le 
calcul des géométries, des fréquences vibrationnelles et autres propriétés des complexes dans les états BS 
et HS, mais elle permet désormais de caractériser leurs énergétiques des états de spin avec une plus grande 
confiance. Des modes opératoires ont été développés à cet effet. Ceux-ci comprennent la reparamétrisation 
des hybrides, l'utilisation de fonctionnelles permettant de descendre l'échelle de Jacob, et la mise en évidence 
d'invariants qui devraient permettre de fortement réduire la dépendance des résultats vis-à-vis du choix de 
la fonctionnelle. Les progrès accomplis concernent également le troisième point de Pauling, et une synergie 
commence à émerger entre les deuxième et troisième points. Ainsi, Les calculs CASPT2 sont maintenant 
en mesure de valider les résultats DFT pour l'énergétique des états de spin des complexes de tailles petites 
et moyennes. A l'avenir, nous nous attendons à un nombre croissant de calculs CASPT2 sur des complexes 
à transition de spin. Toutefois, la DFT et la TDDFT continueront d'être nos outils principaux, notamment 
pour l'étude des systèmes plus grands, et en particulier, la prise en compte des effets d'environnement et 
de photodynamique. Nous avons présenté des premiers travaux allant dans ce sens. 
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